Observations of Water Masers in the NGC 1333 IRAS 4 Region 
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ABSTRACT 



The NGC 1333 IRAS 4 region was observed in the 22 GHz H2O maser hne with an angular 
resolution of about 0.08 arcseconds. Two groups of masers were detected, one near IRAS 4A and 
the other near BI. Among the eight maser spots detected near IRAS 4A, six spots are located 
close to A2, within 100 AU, and the maser velocities are also near the systemic velocity of the 
cloud core. These masers are probably related with the circumstellar disk. Since there is no 
maser spot detected around Al, the star forming process is relatively more active in A2 than in 
Al. Four maser spots were detected near IRAS 4BI. Since most of them are distributed along 
a straight line in the direction of the outflow, BI masers are most likely related with the jet or 
outflow. The disk-outflow dichotomy of II2O masers is discussed briefly. No maser was detected 
near H2 0(C), another maser source reported previously. 



Subject headings: accretion disks 
masers — stars: formation 

INTRODUCTION 



ISM: individual (NGC 1333 IRAS 4) — ISM: jets and outflows 



NGC 1333 is a very active star forming region 
at a distance of 320 pc from the Sun (Jennings 
et al. 1987; de Zeeuw et al. 1999). IRAS 4 is 
a multiple-star system, consisting of at least five 
young stellar objects (YSOs) (Sandell et al. 1991; 
Rodriguez et al. 1999; Looney et al. 2000). IRAS 
4A1 and A2 form a relatively close binary system 
(Lay et al. 1995; Looney et al. 2000). IRAS 
4BI is located ~30" southeast of IRAS 4A. They 
are Class protostars with strong submillimeter 
emission (Sandell et al. 1991). The nature of the 
other objects, IRAS 4BII and IRAS 4C, are less 
clear. 

IRAS 4A is a protobinary system with a sep- 
aration of I'.'S or 580 AU (Looney et al. 2000; 
Reipurth et al. 2002). Al is brighter than A2 in 
the radio continuum emission (Looney et al. 2000; 
Reipurth et al. 2002; Girart et al. 2006; Choi et 
al. 2007). However, the outflow driven by A2 is 
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stronger and larger than that of Al (Choi 2005; 
Choi et al. 2006). The northeastern outflow of A2 
has an intriguing bent feature that was interpreted 
as a result of collision between the outflow and a 
dense cloud core (Choi 2005). Recently, circum- 
stellar disks were discovered by NH3 observations 
(Choi et al. 2007). BI was suspected to be a mul- 
tiple system by Lay et al. (1995). Later, however, 
imaging observations showed BI to be a single ob- 
ject. (Looney et al. 2000; Reipurth et al. 2002). 
The interpretation of Lay et al. (1995) could have 
been a confusion caused by BII. BI has a bipolar 
outflow in the north-south direction (Choi 2001). 

Haschick et al. (1980) detected a water maser 
spot, H20(C), using a single-dish telescope. Since 
H20(C) is located near the bending point of the 
northeastern outflow (see Figure 4 of Choi 2005), 
they could be related. Later, several maser spots 
were detected around A2 and BI (Rodriguez et al. 
2002; Furuya et al. 2003). However, H20(C) was 
not detected by these interferometric observations. 

To study the water masers in the two proto- 
stellar systems, IRAS 4 A and IRAS 4BI, and to 
find the exact location of H20(C), we observed the 
IRAS 4 region using the Very Large Array (VLA) . 
In § 2 we describe our observations. In § 3 we re- 
port the results. In § 4 we discuss the implications 
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of our results and the relation between the water 
masers and the protostar system. A summary is 
given in § 5. 

2. OBSERVATIONS 

The NGC 1333 IRAS 4 region was observed us- 
ing the VLA of the National Radio Astronomy 
ObservatorjQ in the 6ie — > 623 transition of H2O 
(22.235077 MHz). The observations were made in 
two tracks in the A-array configuration. Twenty- 
two antennas were used on 2006 March 24, and 
twenty-three antennas were used on 2006 April 24. 
The spectral windows were set to have 128 chan- 
nels with a channel width of 0.024 MHz, giving a 
velocity resolution of 0.33 km s^^. The phase and 
bandpass calibrator was the quasar 03364-323 (4C 
32.14). The flux calibration was done by observ- 
ing the quasar 0713-^438 (QSO B0710-I-439). The 
flux density of 0713+438 was 0.47 Jy in March 
and 0.44 Jy in April, which are the flux densities 
measured within a day of our observations (VLA 
Calibrator Flux Density Databasell). Comparison 
of the amplitude gave a flux density of 1.00 Jy in 
March and 0.92 Jy in April for 03364-323, which 
agrees with the value in the VLA Calibrator Flux 
Density Database to within 9% and 2%, respec- 
tively. To avoid the degradation of sensitivity ow- 
ing to pointing errors, pointing was referenced by 
observing the calibrators in the V-band (3.6 cm). 
This referenced pointing was performed at the be- 
ginning of each track and just before observing the 
flux calibrator. 

The phase tracking center was a20oo = 03^29™10!529 
and (52000 — 31°13'31'.'000. Imaging was done with 
a uniform weighting, and maps were made using a 
CLEAN algorithm. The synthesized beams were 
FWHM = 0'.'090 X 0'.'075 and P.A. = -41° in 
March, and FWHM = 0'.'094 x 0'.'077 and P.A. 
= —53° in April. The maps have rms noises of 
2.4 and 2.5 mJy beam~^, respectively. The abso- 
lute positional accuracy of the VLA under normal 
conditions is expected to be about O'.'l in the A- 
array configuration. The continuum emission was 
not detected. Figure 1 shows the known compact 
sources in our field of view. 



^The NRAO is a facility of the National Science Founda- 
tion operated under cooperative agreement by Associated 
Universities. Inc. 

^ See |http: / / aips2.nrao.edu/vla/calfiux.html| 



3. RESULTS 

Eight maser spots were detected around IRAS 
4A2, and four spots were detected near IRAS 4BI. 
All the A2 masers are redshifted with respect to 
the systemic velocity of the cloud core (Vlsr = 
6.7 km s-i; Blake et al. 1995; Choi 2001). Out of 
the four maser spots near BI, two spots are red- 
shifted and the others are blueshifted with respect 
to the systemic velocity of the cloud core (Vlsr 
= 6.6 km s^^; Blake et al. 1995). We examined 
the whole field of view, but no maser source was 
detected around H20(C). Table 1 lists the source 
parameters, and Figures 2 and 3 show the distri- 
bution of the detected H2O maser spots. Spectra 
of the H2O masers are shown in Figure 4. 

4. DISCUSSION 

Out of the twelve H2O maser spots detected, 
eight are associated with IRAS 4A2, and four with 
IRAS 4BI. Previous observations with VLA also 
showed several H2O maser spots. BI masers were 
detected by both Rodriguez et al. (2002) and Fu- 
ruya et al. (2003), and A2 masers were detected 
by Furuya et al. (2003) only. 

4.1. IRAS 4A2 

Six maser spots, PC 2-5, 7, and 8, are located 
near the radio continuum source A2, within 100 
AU. The maser velocity ranges from 7.2 to 8.8 
km s~^, and the maser lines are near the systemic 
velocity (6.7 km s~^; Blake et al. 1995; Choi 
2001), within about 2 km s~^. The tight posi- 
tional association with the continuum source and 
the small velocity differences suggest that these 
maser spots may be related with the circumstel- 
lar disk, not the outflow. These maser emission 
probably comes from shocked gas in the accretion 
disk. The exact nature of the shock is not clear. 
There are other examples of H2O masers that are 
probably emitted by the shocked gas in protostel- 
lar disks, including IRAS 00338-^6312, IC 1396N, 
and NGC 2071 IRS 1/3 (Fiebig et al. 1996; Slysh 
et al. 1999; Seth et al. 2002). 

PC 3 and 5 in March are identical to PC 7 and 
8 in April, respectively, in both position and ve- 
locity (Table 1). These two spots seem to have a 
relatively long lifetime, at least a month. Since the 
maser line detected by Furuya et al. (2003) have 
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different velocities, the lifetime is shorter than a 
few years. 

PC 1 and 6 are relatively far from A2, at least 
300 AU. Since they are located in the direction of 
the bipolar outflow, these maser spots are proba- 
bly related with the jet or outflow. 

No maser spot associated with Al was detected, 
which probably indicate that the A2 disk is more 
active than the Al disk. This interpretation is con- 
sistent with the facts that A2 drives more powerful 
outflow than Al and that the A2 disk is brighter 
in NHg lines than the Al disk (Choi 2005; Choi et 
al. 2007). 

4.2. IRAS 4BI 

The systemic velocity of the BI cloud core is 
^6.6 km s""*^ (Blake ct al. 1995), and the maser ve- 
locity ranges from —1.2 to 19.9 km s^^ (Rodriguez 
et al. 2002; this work). The maser spots lie ak)ng a 
O'.'7-long straight line in the southeast-northwest 
direction, P. A. w 29° (Fig. 3). The HCN out- 
flow observed by Choi (2001) flows in a similar 
direction. These facts suggest that the BI masers 
are emitted from the shocked gas related with the 
outflow. Figure 5 shows the position-velocity di- 
agram, which shows that the blueshifted masers 
tend to be located in the southeast direction, while 
the redshifted ones in the northwest. This trend 
agrees with the HCN outflow. 

4.3. Dichotomy between Disk and Outflow 

Torrelles et al. (1997, 1998) claimed that H2O 
masers often prefer to trace selectively either the 
disk or the outflow in a specific region and pro- 
posed that this dichotomy may be caused by the 
differences in the evolutionary stage of the driv- 
ing source. They suggested that the H2O maser 
traces the disk in younger sources and the outflow 
in older sources. 

Our observations appear to show that such a 
dichotomy may exist in the NGC 1333 IRAS 4 
region. Most of the A2 masers trace the disk, 
while the BI masers trace the outflow. However, 
the measured outflow timescales suggest that A2 
(~2000 yr; Choi et al. 2006) is much older than 
BI (~120 yr; Choi 2001). Therefore, the evolu- 
tionary explanation of the dichotomy by Torrelles 
et al. (1997, 1998) does not seem to hold. 

Then what could cause the apparent dichotomy? 



A possible explanation can simply be the low oc- 
currence rate of the maser phenomenon. The H2O 
line amplification requires the existence of warm 
high-density shocked gas (Elitzur 1995), and H2O 
maser in low-mass YSOs is a rare phenomenon. 
Previous surveys of H2O masers toward low-mass 
YSOs showed that the detection rate is typically 
5-20 % (Wilking & Claussen 1987; Terebey et 
al. 1992; Persi et al. 1994; Furuya et al. 2003; 
Gomez ct al. 2006). Since the detection prob- 
ability of either disk-maser or outfiow-maser is 
small, the probability of detecting both disk- and 
outflow-masers around a single YSO would be 
even smaller. For example, let us assume that 
the detection probability of H2O maser from a 
disk is 10 %, and that the probability is similar 
for an outflow. If the two kinds of maser occur 
independently, the probability of detecting both 
in a specific region would be only ^1 %. That is, 
in this example, out of all the regions associated 
with masers, 95 % of the sample would show only 
one kind of masers, either disk-maser or outfiow- 
maser, not both. Therefore, the dichotomy could 
be a result of low probability. 

4.4. H20(C) 

The maser source H20(C) reported by Haschick 
et al. (1980) was not detected in any of the obser- 
vations with VLA (Rodriguez et al. 2002; Furuya 
et al. 2003; this work). There are several possible 
reasons for the nondetection. First, all the VLA 
observations might be unfortunately made dur- 
ing inactive phases of H2 0(C). Second, H2 0(C) 
might be a mis-identification of A2. Haschick et 
al. (1980) observed using a single-dish telescope 
with a large beam size. The right ascension of 
H20(C) coincides with that of A2, and the decli- 
nation of H20(C) could be in error. 

We thank J. Cho for helpful discussions and en- 
couragement. This work was partially supported 
by the LRG Program of KASI. G. P. was partially 
supported by the BK 21 project of the Korean 
Government. 
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Table 1 

H2O Maser Sources in the NGC 1333 IRAS 4 Region 



Source* 










Peak Position 
















AV 


Peak Flux'* 
(Jy beam"*) 




0:2000 




52000 








AS" 






PC 1 


03 


29 


10.190 


31 


13 32.683 


-0.497 ± 0, 


,004 


1 


.683 


± 


0, 


,004 


8.7 


0.7 


0.022 ± 0.002 


PC 2 


03 


29 


10.404 


31 


13 32.454 


-1.598 ± 0, 


,004 


1 


.454 


± 


0, 


,004 


8.8 


0.8 


0.023 ± 0.002 


PC 3 


03 


29 


10.4174 


31 


13 32.3733 


-1.4317 ± 0, 


,0001 


1 


.3733 


± 


0, 


,0001 


8.7 


0.8 


0.735 ± 0.002 


PC 4 


03 


29 


10.401 


31 


13 32.287 


-1.645 ± 0, 


,005 


1 


.287 


± 


0, 


,005 


8.8 


1.1 


0.019 ± 0.002 


PC 5 


03 


29 


10.422 


31 


13 32.246 


-1.373 ± 0, 


,004 


1 


.246 


± 


0, 


,004 


7.3 


0.8 


0.024 ± 0.002 


PC 6 


03 


29 


10.313 


31 


13 30.235 


-2.768 ± 0, 


,005 


-0 


.765 




0, 


,005 


8.8 


1.0 


0.020 ± 0.002 


PC 7 


03 


29 


10.417 


31 


13 32.377 


-1.433 ± 0, 


,001 


1 


.377 


± 


0, 


,001 


8.8 


1.1 


0.121 ± 0.003 


PC 8 


03 


29 


10.422 


31 


13 32.238 


-1.370 ± 0, 


,002 


1 


.238 




0, 


,002 


7.2 


0.8 


0.044 ± 0.003 


PC 9 


03 


29 


12.007 


31 


13 08.008 


18.955 ± 0, 


,004 


-22, 


,992 


± 


0, 


,004 


1.2 


0.7 


0.031 ± 0.003 


PC 10 


03 


29 


12.0128 


31 


13 08.9457 


19.0339 ± 0, 


,0001 


-23, 


,0543 




0, 


,0001 


-1.2 


0.6 


40.973 ± 0.003 


PC 11 


03 


29 


11.9951 


31 


13 08.3301 


18.8073 ± 0, 


,0002 


-22, 


.6699 


± 


0, 


,0002 


19.9 


1.0 


0.588 ± 0.003 


PC 12 


03 


29 


11.937 


31 


13 06.724 


18.057 ± 0, 


,008 


-24, 


.276 


it 


0, 


.008 


19.8 


1.1 


0.024 ± 0.003 



Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and 
arcseconds. All velocities are in km s~*. 

*H20 maser spots presented in this work are labeled with a prefix PC. 

''Position offsets in arcseconds relative to the phase tra<;king center (see § 2). Positional uncertainties were estimated by 
considering the signal-to-noise ratios and the pixel size. 

'^The central velocity and the line width (FWHM) estimated by a Gaussian fit to each spectrum. 

^Fluxes were corrected for the primary beam response. 
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Fig. 1. — Positions of the known compact sources 
within the field of view. Filled circles: Radio con- 
tinuum sources (Looney et al. 2000; Reipurth el 
al. 2002). Plus: Maser source Ha 0(C) (Haschick 
et al. 1980). The size of plus corresponds to the 
uncertainty. 
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Fig. 2. — Map of the H2O maser spots in the 
NGC 1333 IRAS 4A region. Crosses: Maser spots 
detected on March 24. PC source numbers arc 
labeled (see Table 1). The size of markers corre- 
sponds to the synthesized beam: FWHM « O'.'OS. 
Pluses: Maser spots detected on April 24. Open 
circles: Positions of the 3.6 cm continuum sources 
(Reipurth et al. 2002). The straight line at 
the bottom corresponds to 200 AU at a distance 
of 320 pc. Arrows: Direction of the northeast- 
southwestern bipolar outflow of IRAS 4A2 (Choi 
2005). Open ellipse: Schematic diagram of the 
NH3 disk of IRAS 4A2. (Choi et al. 2007). 
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Fig. 3. — Mapof theH20 maser spotsintheNGC 
1333 IRAS 4BI region. Crosses: Mascr spots de- 
tected on March 24. Pluses: Maser spots detected 
on April 24. Filled circles: Maser spots detected 
by Rodriguez et al. (2002). 
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Fig. 4. — Spectra of the H2O masers. (a) Masers 
around IRAS 4A2 from the March 2006 observa- 
tions. (6) Masers around IRAS 4A2 from the April 
2006 observations, (c) Masers around IRAS 4BI 
from the March 2006 observations, (d) Masers 
around IRAS 4BI from the April 2006 observa- 
tions. Vertical dotted line: Systemic velocity of 
the cloud core: Vlsr = 6.7 km s~^ for A2 and 6.6 
km s-i for BI (Blake et al. 1995; Choi 2001). 
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Fig. 5. — Relation between the niascr velocities 
and the angular distances. The horizontal axis is 
the angular separation between the maser spots 
and BI (position of the 3.6 cm source detected 
by Reipurth el al. 2002), projected on a best-fit 
straight line. Pluses: Masers from this work (with 
squares) and from Rodriguez ct al. (2002) (with- 
out squares). The size of the pluses corresponds 
to the uncertainties. Horizontal dotted line: Sys- 
temic velocity of the cloud core. 
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